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Abstract 
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\q , We present the exact formula for neutrino oscillations. By resorting to recent 

results of Quantum Field Theory of fermion mixing, we work out the Green's 
I function formalism for mixed neutrinos. The usual quantum mechanical Pon- 

tecorvo formula is recovered in the relativistic limit. 
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The mixing of mass eigenstates plays a relevant role in high energy physics: the family 
mixing problem of the standard model |TJ], neutrino oscillations and the solar neutrino 
problem || , thermal field theory |4| in hot hadronic matter are some examples where it must 
be considered. However, a thorough theoretical analysis of the mixing of fields in Quantum 
^ \ Field Theory has not been carried out till ref. ||. 

In this paper we present the exact formula for neutrino oscillations: this is done by 
^ ■ studying the Green's functions for mixed fermions. We show that in the study of neutrino 
oscillations it is necessary to take into account the full structure of the Heisenberg vacuum 
for mixed fermions, which exhibits a very rich condensate-coherent state structure |5|.|(| . 

The result is an oscillation formula which differs from the usual one in the non-relativistic 
region, a result which was partially achieved in ref. ||. With the present analysis of the 
full vacuum effects, we get, together with the energy dependent "squeezing" factor of the 
amplitude studied in ref. also an additional term with a different oscillatory frequency: 
resonance is then possible also in vacuum for particular values of the masses and/or of the 
momentum, thus leading to a suppression or to an enhancement of the conversion probability. 

We consider the mixing problem for two Dirac fields with a bilinear interaction La- 
grangian. We denote the fields with v e and (space-time dependence suppressed), referring 
thus explicitly to electron and muon neutrinos: however, the following results are valid for 
general Dirac fields. The Lagrangian we consider has the general form 
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m e ) V e + Vn 



in, 



VeV^ + U U U, 



(1) 



Generalization to a higher number of flavors is straightforward. This Lagrangian can be 
fully diagonalized by substituting for the fields [fj] 



v e {x) = v\{x) cos 9 + z/ 2 (x) sin 6* 
y^{x) = —Pi(x) sin 9 + i>2 (x) cos 9 , 



(2) 



where 9 is the mixing angle and m e = m\ cos 2 9 + m 2 sin 2 9 , = rrii sin 2 9 + m 2 cos 2 9 , 
= (m 2 — rrii) sin 9 cos 9. v \ an d v 2 therefore are non-interacting, free fields, anticom- 



m, 



muting with each other at any space-time point. They are explicitly given by 



e -*>>h,it a r 



1,2 



(3) 



Here and in the following we use t = Xq, when no misunderstanding arises. The vacuum for 
the ai and $ operators is denoted by |0)i j2 : a£ i |0)i j 2 = Pkiltyh 2 = 0. The anticommutation 
relations, completeness and orthonormality relations are the usual ones (see ref. |fl). 

The fields v e and are thus completely determined through eq. @ . In order to circum- 
vent the difficulty of the construction of a Fock space for the mixed fields, which has been 
emphasized in a number of publications, e.g. |7]|§, it is useful to expand the flavor fields v e 
and Un in the same basis as v\ and z/ 2 , 



K(x 



G d \t) v?(x) G e {t) = V- 1 * \<>~ tUlKlt <e(t) + v 



k.v 



,r,a iw hx tp>r\ 



k,l l 



u«{x) = G e l (t)v«{x)G e {t) = V^Y. 



k.v 



k,2 l 



,ikx 



(4) 
(5) 



where Gg(t) = exp 
transformations (|| 

*k,l 



9 f d 3 x ( v\(x)v%{x) — i/2(x)i>i(x) ] is the generator of the mixing 



a. 



Notice that the flavor annihilation operator, say a£ e (£) 
[x) , has contributions also from the anti-particle operator /3 2 ||: 

i£ (t) = cos^i; il + sin^(< 1 < 2 e-^— V M + ^ik^^+^^l,,) . (6) 



The term containing u\^ 1 v s _y L2 -> which is non-zero for mi ^ m 2 , has been wrongly missed in 
the usual treatment of mixing. Relations similar to eq.@ hold for the other operators 0. 

The bilinear mixed term of eq.(|l|) generates four non-zero two point Green's functions 
for the mixed fields u e , v^. The crucial point is about how to compute these propagators: 



if one (naively) uses the vacuum |0)i i2 , one gets an inconsistent result (cf. eq. (|T3D ). Let us 
show this by defining the propagators as 



Sg(x,v) Sf e {x,y) 
Stf(x,y) SfJx,y) 



T[uf{x)vl{y)] T[q(x)V0(y)] 



(7) 



where T denotes time ordering. Use of (fj) gives S ee in momentum representation as 
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# + rni . 2 ft + m 2 



S ee (ko, k) = cos 6 k2 _ ml + iS + sin 6 ~ 2 _ - % - . g , (8) 

which is just the weighted sum of the two propagators for the free fields v\ and v 2 . It coincides 
with the Feynman propagator obtained by resumming (to all orders) the perturbative series 

l-m 2 efl S^S e ) , (9) 

where the "bare" propagators are defined as S e /^ = (ft — m e / ll + id)^ 1 . In a similar way, one 
computes S efM and 5^ e . 

The transition amplitude for an electron neutrino created by a^ e at time t — going 
into the same particle at time t, is given by 

V ee (k, t) = iu^e^ S>(k, f) 7 °<,i , ■ (10) 

where the spinors ui and v\ form the basis in which the field v e is expanded (cf. eq.(^)). Here, 
S> e denotes the unordered Green's function (or Wightman function): iS> e al3 {t, x; 0, y) = 
i,2(0|^e (t, x) z>f (0, y)|0)i j2 . The explicit expression for S> e (k, t) is 

S>^(k, t) = -iJ2 (cos 2 ^ e-M* ^ ojj + sin ^ e -*** ^ ^ ) . (n) 

r 

The amplitude eq. (JT0|) is independent of the spin orientation and given by 

Vee(k,t) = cos 2 9 + sin 2 # |[/ k | 2 e -*K,2-^,i)* . (12) 
For different masses and for k 7^ , |£/k| is always < 1 |J. Explicitly, 
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k,2 "k,! 



|2 



2w M yV 2uJ k,2 J\ (uy + m^uy + mj) 



(13) 



Notice that |£7k| 2 — ► 1 in the relativistic limit |k| ^m\m2 ■ only in this limit the squared 
modulus of P e e(k, t) does reproduce the Pontecorvo oscillation formula. 

Of course, it should be lim t ^ + V ee (t) = 1- Instead, one obtains the unacceptable result 

V ee (k, 0+) = cos 2 # + sin 2 # \U k \ 2 < 1 . (14) 

This means that the choice of the state \0)i t2 i n e q-(0) an d in the computation of the 
Wightman function is not the correct one. The reason is that 1 0} 1^2 is not the vacuum state 
for the flavor fields ||: although the Lagrangian (P is fully diagonalizable by means of the 
rotation (0), this transformation does not leave invariant the vacuum 1 0) 1,2- The mixing 
generator Gg induces on it a SU(2) coherent state structure, resulting in a new state, 

[0(*)>e rf . = G^ 1 (*)IO)i, a , (15) 

which is the flavor vacuum, i.e. the vacuum state for the flavor operators a e //i, fie/n [H : 
Q;^. e , (t)|0(t)) eiAt = /3£ e / M (t)|0(t)) ejA , = 0. It is a condensate of four different fermion- 

antifermion pairs, of the form ct k j/?_ k with i,j = 1,2 (see || for the explicit form). Its 



non-perturbative nature results in the unitary inequivalence with the "perturbative" vac- 
uum lO)^, in the infinite volume limit ||. An important point is represented by the time 
dependence of mixing generator, and consequently of the flavor vacuum. This is not sur- 
prising since the flavor states are not mass eigenstates and therefore the Poincare structure 
of the flavor vacuum is lacking. 

We now show that the correct definition of the Green's function matrix for the fields u e , 
Ufj, is the one which involves the flavor vacuum |0) eiAt , i.e. 

G%(x,y) G${x,y)\ = ( T [^(x)i^(y)] T [v«(x)vg (y)] \ 

Notice that here the time argument y (or, equally well, Xq) of the flavor ground state, is 
chosen to be equal on both sides of the expectation value. Indeed, we observe that transition 
matrix elements of the type ejA1 (0|a e exp [— iHt] c4|0) ejM , where H is the Hamiltonian, do not 
represent physical transition amplitudes: they actually vanish (in the infinite volume limit) 
due to the unitary inequivalence of flavor vacua at different times (see below). Therefore 
the comparison of states at different times necessitates a parallel transport of these states 
to a common point of reference. The definition ( |16|) includes this concept of parallel trans- 
port, which is a sort of "gauge fixing": a rich geometric structure underlying the mixing 
transformations (|j) is thus uncovered. In a forthcoming publication we will discuss these 
geometrical aspects of field mixing, including Berry phase and gauge structure. 
In the case of v e — > v e propagation, we now have (for k = (0, 0, |k|)): 

G ee (k ,k) = S ee (k ,k) + 2n i sin 2 9\\V k \ 2 (ft + m 2 ) 5(k 2 -m 2 2 ) 

- \Uk\\V k \ ^ ^-k,2 s (ka ~ w 2 ) + e r ^- k ,2 K,2 $( k o + ^2)) , (17) 



where we used e r = (— l) r and |Vk| = yl — \Uk\ 2 . Comparison of eq.(pT7|) with eq.® shows 
that the difference between the full and the perturbative propagators is in the imaginary 
part. 

The Wightman functions for an electron neutrino are iG^^(t, x; 0, y) = 
e , M (0|i£(*,x) ^(0,y)|0)e lA ,, and iG>f(t, x; 0, y) = e , M (0|z/«(f, x) pf(0, y)|0) e>M . These are 
conveniently expressed in terms of anticommutators at different times as 

^(k, t) = J2 <i {<M <} ^ l + u -k,i <} ^} , (is) 

r 

*G>f%t) = £ [<Ki {« k „(t),<}e-^ t + ^, 2 < {P r \^)><}e iu>k > 2t ] ■ (19) 

r 

Here and in the following a^ e stands for a^ e (0). We now have four distinct transition 
amplitudes, given by anticommutators of flavor operators at different times: 

V r ee (k, t) = 1 ti&e*-"* G> (k, t) 7 X,i = {<, e (t), <} 

= cos 2 6> + sin 2 6> [|[/ k |2 e -iK 2 -^,i)* + |y k |V( w M+"*,i)*] , (20) 
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= e r |^k||^ k | sin 2 fl [e^^,^ _ e -<(«*, 2 +«fc,i)*] , (21) 

V; e (k,t) = i«* e **.'* G > (k,t) 7 V k)1 = { <,,(*), <} 

= |C/ k | cos# sin6> [l - e *( w fc,a- w *,i)*] , (22) 

7^(k,t) EE ^V-^G> (k,t) 7 V kjl = {CUW,< e } 

= e r |Vk] cos# sin# [l - e'^^^] . (23) 

All other anticommutators with vanish. The probability amplitude is now correctly 
normalized: lim t ->o+Vee(k,t) = 1, and P ge , T^e, ^//e go to zero in the same limit t — > + . 
Moreover, 

in r e (k,t)| 2 + |PL(k,t)| 2 +|^ e (k,t)| 2 +|pr ( k ,t)| 2 = l, (24) 



as the conservation of the total probability requires. We also note that the above transition 
probabilities are independent of the spin orientation. 

For notational simplicity, we now drop the momentum and spin indices. The momentum 
is taken to be aligned along the quantization axis, k = (0, 0, |k|). It is also understood that 
antiparticles carry opposite momentum to that of the particles. At time t — 0, the vacuum 
state is |0) 6i/1 and the one electron neutrino state is 

\u e ) = al\0) e>lt = cos0a\ + \U\ sin# a\ - e \V\ sin# a\a\p\ |0> 1;2 - (25) 

In this state a multiparticle component is present, disappearing in the relativistic limit 
|k| ^> y/mim2 '■ in this limit the (quantum-mechanical) Pontecorvo state is recovered. The 
time evoluted of \v e ) is given by \v e {t)) = e~ lHt \i , e ). 

Notice however that the flavor vacuum \0) e ^ is not eigenstate of the free Hamilto- 
nian H. It "rotates" under the action of the time evolution generator: one indeed finds 
limy^oo e ,^(0 | 0(t)) etfM = 0. Thus at different times we have unitarily inequivalent flavor 
vacua (in the limit V — > oo): this expresses the different particle content of these (coherent) 
states and it is direct consequence of the fact that flavor states are not mass eigenstates. 

As already observed, this implies that we cannot directly compare flavor states at different 
times. However we can consider the flavor charge operators, defined as Q e /^ = a \/^ a e/n — 

Pl/pPe/p- We tneI1 haVe 

w(0(t)\QMt))e# = e,Mt)\QMt))e,» = 0, (26) 
(Ve(t)\QeMt)) = | {«.(*) , 4 } \ ' + | (t) , a\ } | 2 , (27) 

(v e {t)\Q»\v e {t)) = |K(t),«l}| 2 + M(t),al}\ 2 . (28) 

Charge conservation is ensured at any time: (f e (i)| (Qe + Qp) \ u e{t)) = 1 and the 
oscillation formula readily follows as 
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(29) 



= 1 - sin 2 (26>) 



\Uvr sin' 



t + Vk sin' 



2 • 2 / ^,2 + ^,1 



= sin 2 (26>) 



(30) 



2 • 2 / W fe,2 — ^fe,l 



t4 sin 2 



i + iVLr sin 2 



2 / w fc,2 + ^fe,l 



This result is exact. There are two differences with respect to the usual formula for 
neutrino oscillations: the amplitudes are energy dependent, and there is an additional os- 
cillating term. For |k| ^> yjm\m<i-, |^k| 2 — ► 1 and |Vk| 2 — > and the traditional oscillation 
formula is recovered. However, also in this case we have that the neutrino state remains 
a coherent state, thus phenomenological implications of our analysis are possible also for 
relativistic neutrinos. Further work in this direction is in progress. 

This work is partially supported by MURST, INFN, and by a network supported by the 
European Science Foundation. 
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